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We study systematically an extended Bose-Hubbard model on the triangular lattice by means of 
a meanfield method based on the Gutzwiller ansatz. Pair hopping terms are explicitly included and 
a three-body constraint is applied. The zero-temperature phase diagram and a variety of quantum 
phase transitions are investigated in great detail. In particular, we show the existence and stability 
of the pair supersolid phase. 
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I. INTRODUCTION 

The development in experimentally manipulating 
ultra-cold atoms in an optical lattice has allowed the 
realization of novel quantum states and quantum phase 
transition in strongly correlated systems e.g., a quan- 
tum phase transition from a superfluid (SF) to a Mott 
insulator (MI) has been predicted and observed Q. The 
condensation of paired electrons, which provides the basis 
of superconductivity of metallic superconductor, plays an 
essential role in modern condensed-matter physics. Thus 
realizing pairing related novel quantum states in the con- 
text of ultra-cold atoms has attracted considerable recent 
interest, both in theoretical and experimental research. 

It has been demonstrated recently that such states can 
be realized for lattice bosons with attractive on-site in- 
teractions, which is stabilized by a three-body constraint 
0, The three-body constraint has been realized by 
large three-body loss processes d, H|. The system can 
be mapped into spin-1 atoms at unit filling [7]. Besides 
the conventional atomic superfluid (ASF) phase, a pair 
(dimer) superfluid (PSF) phase consisting of the conden- 
sation of boson pairs emerges under sufficiently strong 
attraction^, 4]. The PSF state is manifested as a second- 
order processes of the single-atom hopping in the optical 
lattice. Various phase transitions among the ASF, MI 
and PSF are investigated in great detail QMIsHIIi]- Both 
the ground-state and the thermal phase diagrams are ob- 
tained. However, the single-species pair supersolid (PSS) 
was not found in the present system when the nearest- 
neighbor repulsion is included, except for an isolated con- 
tinuous supersolid at the Dirac point |3] . The three-body 
constraint and the effective pair hopping makes the sys- 
tem resemble hardcore bosons when onsite repulsion U is 
weak. Thus the PSS state might suffer from the same in- 
stability of the supersolid state (SS) for hardcore bosons 
on the square lattice [l3|, [T3] . The pair supersolid state 
was predicted only in system with correlated hopping 
[TH and in bilayer system [TtJ ■ 



For hardcore bosons on the triangular lattice, su- 
persolid state emerges basing on an order-by-disorder 
mechanism , by which a quantum system avoids classical 
frustration |18M2l| . Aiming to realize PSS state, we thus 
focus on lattice bosons on the triangular optical lattice. 
The system we consider is an extended Bose-Hubbard 
model with the three-body constraint a] 3 = 0. The pair 
hopping terms are explicitly included. The Hamiltonian 
is 

H = -t^2(a\aj + OjOi) - t p ^2 l {a\a\ajaj + aJ^aidi) 



+ ^-^2ni(ni - 1) + V'^n i nj - fi^n,., 



(1) 
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where denotes nearest-neighbor sites, a\{a.i) is the 
boson creation (annihilation) operator at site i, and n, = 
ajflj the boson number operator; t (t p ) is the single-atom 
(pair) hopping amplitude, U the on-site interaction, /i the 
chemical potential, and V the nearest-neighbor repulsion. 
Such a system can be realized experimentally for dipolar 
bosons polarized by an external electric field and confined 
in an optical lattice (22j . The pure pair hopping limit (t = 
0) can be realized by a mechanism based on collisions that 
induce transport [23J. At the limit V = 0, the system 
can be realized in an atom-molecule coupling system in 
a state-dependent optical lattice 24] . A similar effective 
model of bosons with the three-body constraint can be 
realized with spin-1 atoms 

In present work we systematically study the model by 
means of the meanfield method based on the Gutzwiller 
ansatz (25l - [28j . The zero-temperature phase diagram of 
the system is studied in great detail. The existence of the 
PSS phase at the t — > limit is explained by a mapping 
between the present model and the hardcore bosons on 
the triangular lattice. We then show that the PSS phase 
is stable in the presence of small on-site repulsion U and 
single-atom hopping t. 

This paper is organized as follows: we first discuss the 
classical limit of model ([T]) at zero temperature in SeclHl 
The model shows various solid states, which are the ba- 
sis to form a pair supersolid state. Then we describe the 
mean-field method in Sec. IIII1 We present our main re- 
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suits in Sec. [IVJ We discuss the the noninteracting case 
(U = V = 0) in Sec. IIV Al For interacting cases, we 
present the zero-temperature phase diagram in the limit 
single- atom hopping t — 0, focusing on the parameter 
region where the pair supersolid phase emerges, in Sec. 
IIVBI We then investigate if the results are stable against 
putting a finite t/t p in Sec. IIV C[ in which a more re- 
alistic phase diagram is presented. We conclude in Sec. 

13 



II. CLASSICAL LIMIT 

There are three sublattices A, B and C in the triangu- 
lar lattice, as shown in Fig. [T] (a). In the classical limit 
(i = 0, t p = 0) and zero temperature, the energy per unit 
cell is 

£ 'a= ^{n A {nA - 1) + n B {n B - 1) + n c (n c - 1)) 
zV 

- K n A + n B + n c ) + — (n A n B + n B n c + n c n A J,2) 

where z = 6 is the coordination number of the triangu- 
lar lattice. nA,n B ,nc is the occupation number on the 
sublattice i = A, B, and C, respectively. 

A solid state is formed when the symmetry of occu- 
pation on the three sublattices is broken spontaneously. 
The order can be represented by the occupation pattern 
(n A ,n B ,nc), with n A ,n B ,nc not equal. If n A — n B = 
nc, the pattern represents an MI state. 

By comparing the energy per unit cell, we obtain the 
phase boundaries between various ordered phases, as 
shown in Fig. [T^b). The solid states (0,0,2), (0,1,2), 
(0, 2, 2), (1, 1, 2) and (1, 2, 2), which have two bosons sit- 
ting on one or two sublattices, are of special interest. 
The corresponding density is p = 2/3,1,4/3,4/3,5/3, 
respectively. Other solid states related with patterns 
(0, 0, 1), (0, 1, 1) and MI states, (0, 0, 0), (1, 1, 1), (2, 2, 2), 
with corresponding densities p = 0, 1, 2, are also found. 




(a) (b) 

FIG. 1. (Color online) (a) The triangular lattice and its 
three sublattices. The (0, 1, 2) solid is also shown. (b)Zero- 
temperature phase diagram in the classical limit t — 0, t p — 0. 



III. MEANFIELD METHOD 

We start with the Gutzwiller ansatz of the many-body 
wave function 

i^>=n(i>^' n >)' ( 3 ) 

i n 

where \i,n) is the Fock state of n bosons occupying the 
site i, and Ci tTl is the probability amplitude. The three- 
body constraint requires the maximum number of bosons 
on each site n max = 2. The mean- field decoupling ap- 
proximation 

a i a J = ( a i) a J + a l( a j) ~ ( a l)( a j) , ( 4 ) 

= (4 2 ) a j + a l 2 ( a j) ~ ( a i 2 )( a |) > ( 5 ) 

mrij = (rii)nj + ni{rij) - (m)(nj) (6) 

are exact for such a state. The atomic superfluid order 
parameter is defined as 4> a ,i = (a>i) = c *o c i,i + V^c^Cj^ 
and the pair superfluid order parameter is 4> Pi i = (af) = 

V2c* a : 2- The average occupancy (m) = J2n n \ c i:-n\ 2 ■ 
For simplicity, we choose the order parameters to be real, 
i.e., 4> a (p).i — ^a(p) v The meanfield version of the Hamil- 
tonian ([1]) is thus written as a sum over single-site terms 
H = 2^, with 

Hi = -t{4> a ,ia\ + H.c.) - t p (4> p4 al 2 + H.c.) 

+ynj(rii - 1) + V(nifii - n \ - p^ (7) 

where <f> a ( P )A = E(,),^( P )j and = Eo>,K) sum 
over sites j neighboring to site i. Written in the matrix 
form, 

(d(0) -t<t> a ,i -V2tp<j> p>i \ 

-t4> a< i d(l) -\/2#a,i , (8) 

-V2tp4>p t i -V2hf> a> i d(2) ) 

where the diagonal element 

d(k) = ^k(k - 1) + V(k - K -^)n t - fik 

(9) 

with k = 0,1,2. 

Without the nearest-neighbor repulsion (V = 0), there 
is no solid ordering. The Gutzwiller ansatz §5§ reduces 
to 

|^) = ( Co |0)+ Cl |l)+ C2 |2)) Ar , (10) 

with N the number of lattice sites. Therefore, 4>a(p),i — 

<f>a(p),( n i) = ( n )> and 4>a(p),i = z<j) a ( p ),n t = z(n). We 
then find the ground state in a self-consistent way: Given 
an initial state |i/>), the order parameters <p a r p ) and the 
occupancy (n) are calculated. The matrix (|5J) is thus ob- 
tained and diagonalized. 4> a ( p ) and (n) are then evaluated 
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again from the ground state and put back to ©. This 
is done iteratively until the order parameters converge to 
a self-consistent solution. There can be more than one 
self-consistent solutions. The one with lowest energy is 
chosen. This procedure is equivalent to minimization of 
the ground energy with respect to the order parameters 
25]. 

With the nearest-neighbor repulsion (V ^ 0) turning 
on, a solid order may emerge. To accommodate the solid 
order, we assume 



IV;) = 11 (c i ,o|»,0) + c i)1 |»,l) + c i)2 |t,2))"S (11) 

i=A,B,C 

where Ni is the number of sites in sublattice i £ A, B, C. 
The solid order is represented by the occupation pattern 
((ua), (tib), ( n c))i which can be depicted by a solid order 
parameter Ap 2 = J2(( n i) ~ P) 2 > wrtn p — |XX n ») tne 

i i 

density of bosons. For a site in sublattice A, one has 



<Pa(p),A = ^(4'a( P ),B + 4>a(p),c), 

Z 

™a = tjOb) + (nc))- 



(12) 



Their equivalences under cyclic transformation can also 
be found. The ground state can be found self- 
consistently. We start with an arbitrary state, from 
which the parameters (j) a ( p ) f A, <Pa(p),B, <t>aM,c> ( u a), 
(tib), and (nc) are evaluated. The matrix (J8j for sub- 
lattice A is first obtained and diagonalized. From the 
obtained ground state, "new" <fi a ( p ),A and (tia) are cal- 
culated. Then the matrix © for sublattice B is obtained 
and diagonalized. "new" 4> a { P ),B and (tib) are obtained 
and used in the matrix © for sublattice C. This is done 
recursively, until the estimated order parameters and oc- 
cupancies converge. Again, the lowest energy solution is 
chosen. 

The zero-temperature phase diagram of the system is 
constructed by checking the obtained solutions. The MI 
and solid states are characterized by <p a = <fi P = 0, 
and Ap 2 = and Ap 2 ^ 0, respectively. The ASF 
and atomic super-solid (ASS) phases are characterized 
by (f> a ^ 0, 4> v ^ 0, and Ap 2 = and Ap 2 ^ 0, respec- 
tively The PSF and PSS both have the pair superfluid 
character <f> a = and <j) p ^ 0. Ap 2 = in the PSF 
state, but Ap 2 ^ in the PSS phase. In the PSS state, 
<Pp,a, 4>p,b and tfi Pi c arc different in general. We take the 
mean value <f> p = J2i 4>p,i/3 as the order parameter. 



IV. RESULTS 



In this section, we present our main results. 



A. Competition between t and t p in the 
non-interacting case 



We start with the non- interacting limit, focusing on 
the competition between the single-atom hopping t and 
the pair hopping t p . Without U and V, the diagonal 



element ^ of ([5} further reduces to 

d(k) = —pk + zt<f) 2 a + zt p (f) 2 p . 



(13) 



The meanfield solutions for various p/t and t p /t are found 
in the way described in Sec. IIIII for V = case. The 
phase diagram is constructed by checking the obtained 
solutions, as shown in Fig. [2] (a). There is a p = MI 
(empty) phase and a p = 2 MI phase at negative and 
large chemical potentials, respectively The MI phase 
is the result of the three-body constraint. At small pair 
hopping, there is an ASF phase between the empty phase 
and the MI phase. The transitions between the ASF 
phase and the two MI phases can be continuous or first 
order, depending on the ratio t p /t. This is because, in 
presence of t p , an energy level crossing may preempt the 
continuous evolving of the ground state. The first or- 
der transition behavior is shown in Fig. [5] (b), in which 
two energy levels cross each other, leading to a discontin- 
uous change of order parameters. Similar MI-ASF transi- 
tion behaviors are found in the attractive Bose-Hubbard 
model with three-body constraint Q . As pair hopping t p 
becomes large enough, an ASF-PSF transition happens 
as expected. The transition is found to be first order, or 
continuous, depending t p /t, as illustrated in Fig. [5] (a). 
Several tricritical points where the jump of order param- 
eters vanishes are found. However, according to [3J, |9(, 
the coupling of the Ising order parameter and the Gold- 
stone mode will drive this quantum phase transition to be 
a weakly first-order one through the Coleman- Weinberg 
mechanism [29j . For such a weak first-order transition 
the mean- field theory may be not reliable. 
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FIG. 2. (Color online) (a) Zero-temperature phase diagram 
at U — 0, V = 0. The dashed (blue) lines and the black 
dot represent continuous phase transitions. The solid (red) 
lines represent first order transitions. The open circles denote 
the tricritical points. The insets show the details around the 
ASF-PSF tricritical points, (b) Levelcrossing and the ASF 
order parameter <f) a jump at the MI (p — 0)-ASF transition 
point, along the t p /t = 1 line. 
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B. The limit single-atom hopping t = 

We now turn to the interacting case. To demonstrate 
the effect of the pair hopping term, we first discuss the 
limit t = 0, at which Eq. (JSJ reduces to 



d(0) -V2t p (f> p ,, 

II; ! (2(1) 

-V2t P 4> P!i d(2) 
where the diagonal elements read 



d(k) = ^k(k-l) + V(k-^)fhi- pLk+t p <f> p , i( j> P}i , (15) 

with k — 0,1,2. Without t, the Hamiltonian bears the 
particle-hole symmetry, which means that it is invariant 
under the transformation n — > 2—n and p — > 2zV—p+U. 



1. t = 0, U=£Q, V = 

With the on-site repulsion U present, the PSF state is 
the only possible superfluid state and a p = 1 MI state 
emerges. Without the nearest-neighbor interaction V, no 
solid order presents. <p p> i — 4> p is uniform and (j) Pi — z<fi p . 
Competition between the on-site repulsion and the pair 
hopping leads to phase transitions among MI and PSF 
phases. 

The diagonal terms in Eq. ([14]) further reduces to 
d(k) — ¥^k(k — 1) — pk + zt p (/) p , and the three eigenvalues 
are found as follows 



zt^l + ^-n : \ {- -//)-' + 2( -.1 n )in f,. 



U 
~2 



E 3 = zt p (f) p - p 



(16) 



Ei is the ground state if p < or p > U, consider- 
ing U > 0. Note that (j> p is self-consistently determined 
by the eigenvector v\ = (co,0, c-i) associated to E\ as 
4> p = V%C2Co and the eigenvector associated to E 3 is 
V3 = (0,1,0). Given parameters t p ,p and U, the self- 
consistent solution E\{(j) p ) should satisfy the condition 
dEi/d<f> p — 0, which determines the superfluid parame- 
ter 4> p of the ground state: 



p)/zt p 



'IL 

(A 



((f - M )Atp) 2 >i; 



H)l 'ztpf < 1. 



(17) 

Thus the PSF and the MI (p = 2) boundary is p/U - 
1/2 = ztp/U, where Vi evolves to (0, 0, 1) and </> p = 0. 
The PSF and the MI (p = 0) boundary is p/U -1/2 = 
—ztp/U, where V\ evolves to (1, 0, 0) and <j> p = 0. 

If < p/U < 1, then the eigenlevels E\ and E3 may 
cross each other at a first order transition point. If both 
sides of the point has (j> p = 0, then E\ = E3 yields p/U — 
1, which is the MI (p = 1) and MI (p = 2) transition line, 
or, p/U = 0, which is the MI (p = 1) and MI (p = 0) 
transition line. It is still possible that 4> p ^ on one side 



of the transition and (j) p = on the other side, which 
leads to the following equation 



^p©p ~ I - /' - \/(| - J /'- l|Si 



with 4> p ^ given by Eq. ([Tf). This gives the PSF and 



(14) MI (p 



1) transition line 
2zt p /U = 1 + y/T 



(2p/U - l)^ 



(19) 



The phase diagram is shown in Fig. [3Ja), which can 
be compared with the similar phase diagram presented 
in Ref. [3], where the three-body constraint was not 
applied. Due to the particle-hole symmetry, the phase 
diagram is symmetric about p/U = 1/2. The nature 
of phase transitions involved is further investigated, as 
illustrated in Fig. Elb), where p/U varies from -1 to 2 
along the line t p /U = 0.15. It is clear that the transitions 
between the p = 0, p = 2 MI phases and the PSF phase 
are continuous, while the transition between the p = 1 
MI and the PSF is first order. The latter behavior is 
different from the nature of the SF-MI transition in the 
softcore Bose-Hubbard model with single-atom hopping. 
The p = 1 MI state v 3 = (0, 1,0) can not evolve to the 
PSF state v\ continuously in the truncated Hilbert space 
due to the three-body constraint. 
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FIG. 3. (Color online) (a) Zero-temperature phase diagram 
at t = 0, V = 0. The solid (red) lines denote first-order phase 
transitions, while the dashed (blue) lines represent continuous 
phase transitions, (b) Density p and pair superfluid parameter 
4> v versus p/U at t p /U — 0.15. 



2. t = 0, U = 0, V 

With the nearest neighbor interaction V turning on, 
but the on-site repulsion turning off, two solid phases: (0, 
2, 2) and (2, 0, 0), emerge, as shown in Figf4j Between 
the two solids, there are two PSS phases (A and B) which 
are characterized by local density fluctuations Ap 7^ 
and the mean PSF order parameter 4> v ^ and the ASF 
order parameter <f) a — 0. The boson density in PSS-A p < 
1, while that in PSS-B p > 1. The PSF phase exists at 
the outside of the two solids and the PSS phases, between 
the empty phase and the p = 2 MI phase. Due to the 
particle-hole symmetry, the phase diagram is symmetric 
about p/V = 6. 
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The phase diagram is much like that of the meanficld 
hardcore Bose-Hubbard model on the triangular lattice 
30] . The only difference is that the pair states replace the 
corresponding single-atom states: the PSS state replaces 
the ASS state and the PSF state replaces the ASF state. 
This can be understood in the following way. The matrix 
of the meanfield Hamiltonian of the hardcore bosons on 
the triangular lattice is 



H 



(he) 



d(0) 



~t4>a,; 
d(l) 



where the diagonal element is 



d(k) = V(k 



-)rii - flk + t$a,i§a, : i 



(20) 



(21) 



with k = 0, 1. Considering [7 = 0, the ground state of Hi 
in Eq. (TT4"]) is always in the subspace extended by the two 
Fork states |0) and |2). The matrix of Hi is thus equiva- 
lent to Hf 1 ^ after applying the mapping 2t p — > t, 2\x — > 
li. The ground state of H^. \ip)i = c i |0) 4- Cj 2 1 2) maps 

to the ground state of fff lc) : |^ (hc) ) 2 = c i)0 |0) + c 2 , 2 |l), 
which leads to <f> p — \/2<p a , with (f> p the PSF order param- 
eter of Hi and <f> a the ASF order parameter of hardcore 
bosons. Therefore, the phase diagram Fig. [3] can be ob- 
tained from the meanfield phase diagram of the hardcore 
bosons on the triangular lattice |30| after applying the 
mapping. 
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FIG. 4. (Color online) Zero-temperature phase diagram at 
t = 0, U — 0. The solid (red) lines denote the first-order 
phase transitions, while the dashed (blue) lines and the solid 
(blue) dot represent continuous phase transitions. 

The properties of phase transitions are also obtained. 
Transitions from the PSF to the two solids((0,0,2) and 
(0,2,2)), and to the PSS states are first order, but those 
from the PSF to the MI phases and from the PSS to the 
two solids are continuous, except for the \i/V — 6 point 
marked by a blue dot in Fig|4] [H], [32| . The transition 
from PSS-A to PSS-B is first-order. These results are 
consistent with what found for hardcore bosons on the 
triangular lattice [Til HH, In a recent work[33j], a 
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FIG. 5. (Color online) (a),(c),(d) Zero-temperature phase di- 
agram for U/V — 2,3.1,4, respectively. First-order phase 
transitions are denoted by solid (red) lines, whereas the con- 
tinuous phase transitions are represented by dashed (blue) 
lines and the solid (blue) dot. (b) The density p, PSF order 
parameter <f> p , and solid order parameter Ap versus p/V at 
t p /V = 0.36, U/V = 2. 



Quantum Monte Carlo study also predicts the PSS phase 
between the two solids. However, the PSS region is much 
smaller due to quantum fluctuations. 



3. t = 0, U # 0, V 

With the on-site repulsion U present, more phases 
emerge, as shown in Fig. [5] When U/V is not large 
(< 3.3), the PSS phase persists. In the strength U/V > 
3, a p = 1 MI phase (1, 1, 1) emerges and the solid 
phases (1, 1, 2) and (0, 1, 1) take over the (0, 2, 2) and 
(0, 0, 2) solids. Nevertheless, as long as U/V < 3.3, the 
pair hopping term can still lower the ground state en- 
ergy by forming the PSS-A or B states on the basis of 
(2, 0, 0) or (0, 2, 2) solid, respectively. For example, at 
U/V = 3.1, fi/V = 6, t p /V = 0.48, we find the ground 
state is a PSS-A state 

\tp) = (0.353|0) + 0.93512))^ (0.912|0) + 0.40912))^ 
(0.912|0) + 0.40912))^°, (22) 

which leads to the occupancy (1.750, 0.335, 0.335), ASF 
order parameter <f> a = 0, and the mean PSF order pa- 
rameter <fi p = 0.508. Above the threshold U/V = 3.3, 
the PSS phases disappear finally. The properties of asso- 
ciated phase transitions is investigated by checking vari- 
ous order parameters, as illustrated in Fig. [Sfb), in which 
various order parameters are plotted as functions of pi/V 
at U/V = 2 1 t p /V = 0.36. 
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C. Finite t/t p in the interacting cases 

We now turn to discuss if the above results are stable 
against putting in the single-atom hopping t. 



1. V = 0,U=£0 

We first study the change of the phase diagram shown 
in Fig. El 

For small enough t, two ASF phases emerge at the 
boundary of the p — 1 MI and the p = 2 MI phases, and 
at that of the p = 1 MI and the p — MI phases, as shown 
in Fig. [Ha). The MI (p = 2)-PSF and the MI (p = 0)- 
PSF transitions are continuous, while the MI (p = 1)- 
PSF transition is still first order. The MI-ASF transi- 
tions and the ASF-PSF transition changes from continu- 
ous to first order, due to the mechanism that levelcorss- 
ing preempts the continuous evolving, according to our 
meanfield analysis. Here the aforementioned coupling of 
the Ising order parameter and the Goldstone mode may 
change phase transition behavior. Further investigation 
beyond the mean-field theory is needed. 



1.5 

I ' 

0.5 


-0.5 



(a) MI(p=2) 

— t/V =0.02 
ASF 



PSF 




0.04- 



ASF 



i PSF 



MI(p=0f 



0.2 



tp/U 



FIG. 6. (Color online) Zero-temperature phase diagram for 
V = at t/U = 0.02 (a) and at t/U = 0.04 (b). The 
solid (red) line denotes first-order phase transition, while the 
dashed (blue) lines represent continuous phase transitions. 
Open circles denote tricritical points. 

When t is large enough, the p = 1 MI phase is ex- 
cluded by the ASF phase completely, as shown in Fig. 
|6] (b). Quantum hopping dominates the physics. The 
phase diagram shows the same topology as that for non- 
interacting case (Fig. [2ja)). 
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FIG. 7. Zero-temperature phase diagram at U = 0, t/V = 
0.15 (a), t/V = 0.25 (c) and t/V = 0.4 (d). The solid (red) 
lines denote the first-order phase transitions, while the dashed 
(blue) lines and the solid (blue) point represent continuous 
phase transitions. The open circles denotes tricritical points. 
Levelcrossing and the jumps of order parameters of the PSS-A 
to PSF transition at p/V — 4.5 is shown in (b). 



The ASS-PSS transition is continuous and the PSS-PSF 
transition is first order, as illustrated in Fig. [7]Jb). 

When t becomes stronger, the ASF region is enlarged 
and encloses the solid and the ASS phases, and the PSS 
phases disappear, as shown in Fig. [7{c). One interest- 
ing phenomenon is that the ASS-A to the (0, 0, 2) solid 
phase changes from first order to continuous when t p /V 
is weaken, while the ASS-B to the (0, 2, 2) solid phase is 
always continuous. 

When t becomes even stronger, all solid phases and 
ASS phases disappear. The dominate physics is the ASF- 
PSF phase transition. The phase diagram for t/V = 0.4 
is presented in Fig. EJd), showing the same topology as 
the non-interacting case (Fig. [UJa)). 



3. Both U / and V 



2. U = butV 

We now check the stability of the phase diagram Fig.0] 
in the presence of the single-atom hopping t. 

With small t present, an ASF region emerges between 
the p — 2 MI and the (0, 2, 2) solid phases, meanwhile 
another ASF phase emerges between the p = MI and 
the (2, 0, 0) solid phases, as shown in Fig. [7{a). In 
addition, two ASS phases (A: p < 1, B: p > 1) emerge in 
the PSS region at small t p /V, between two solid phases. 
Large PSS region persists in the region with larger t p /V. 



In a real experimental system, the quantum hoppings 
and atom interactions can exist simutaneously. The in- 
terplay among them results in complex, but interesting 
physics effects. We fix the ratios pt/V = 5.8 and t p /t = 5 
and show the zero-temperature phase diagram in the pa- 
rameter space tp/V and U/V, see Fig. |5{a). 

When the nearest-neighbor repulsion V dominants, the 
system is in the (0, 0, 2) solid phase. As hopping strength 
grows, the system goes to the PSS-A state through a con- 
tinuous phase transition. Further enlarging the hopping 
strength results in the PSF state through a first order 
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transition. The PSS-A phase is stabilized in a large pa- 
rameter region. 

On the other hand, the system is in the (0, 1, 1) solid 
state, if the on-site repulsion U dominants. As hopping 
strength increases, the system undergoes a continuous 
transition to the ASS-B phase. By increasing further the 
hopping strength, the solid order is destroyed and bosons 
undergo a first order transition to the ASF phase. Even 
larger hopping with the fixed ratio t p /t — 5 drives the 
system into the PSF phase through a continuous transi- 
tion, as expected. 

In the middle value of U /V and t p /V, a small region 
of ASS-A phase exists. 



1 ' \ , 
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FIG. 8. (Color online) Zero-temperature phase diagram at 
tp/t = 5,fi/V = 5.8. The solid (red) lines denote the first- 
order phase transitions, while the dashed (blue) lines repre- 
sent continuous phase transitions. 



V. CONCLUSION 

We have studied systematically the extended Bose- 
Hubbard model on the triangular lattice, in which the 
pair hopping terms are explicitly included and the three- 
body constraint is applied, by means of mean-field ap- 
proaches based on the Gutzwiller ansatz. The zero- 
temperature phase diagram and various quantum phase 
transitions are investigated in great detail. In particular, 
the existence and stability of the pair supersolid phases 
are shown. At the limit that the single- atom hopping 
is zero, we provided the mapping between the present 
model and the hardcore bosons on the triangular lattice. 
The existence of the PSS phase is thus understood. We 
have also shown that the PSS phase are stable under 
the perturbation of the on-site repulsion and the finite 
single-atom hopping. Experimentally, the triangular op- 
tical lattice can be implemented [34| . The three-body 
constraint can be realized with spin-1 atoms or large 
three-body loss processes [ij @. Our results for the pure 
pair hopping limit (t — 0) are applicable to the system 
in which bosons pair hopping based on collisions that in- 
duce transport [23|. The results in the parameter space 
V = are applicable to the atom-molecule coupling sys- 
tem proposed in (24|, when the three-body constaint is 
applied. For V 7^ 0, our results are useful to analyse 
dipolar bosons polarized by an external electric field and 
confined in an optical lattice [22| ■ However, in such a 
system, the occupation-dependent single-atom hopping 
needs to be discussed. 
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